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An example of computation of the density of 
ones in probabilistic cellular automata by direct 
recursion 


Henryk Fuks 


Abstract We present a method for computing probability of occurence of Is in 
a configuration obtained by iteration of a probabilistic cellular automata (PCA), 
starting from a random initial configuration. If the PCA is sufficiently simple, one 
can construct a set of words (or blocks of symbols) which is complete, meaning 
that probabilities of occurence of words from this set can be expressed as linear 
combinations of probabilities of occurence of these words at the previous time step. 
One can then setup and solve a recursion for block probabilities. We demonstrate 
an example of such PCA, which can be viewed as a simple model of diffusion of 
information or spread of rumors. Expressions for the density of ones are obtained 
for this rule using the proposed method. 


1 Introduction 

Binary probabilistic cellular automata (PCA) in one dimension are one of the most 
frequently studied types of cellular automata, and one of the most natural and most 
frequently encountered problems in PCA is what the author proposes to call the 
density response problem: If the proportion of ones in the initial configuration drawn 
from a Bernoulli distribution is po, what is the expected proportion of ones after t 
iterations of the PCA rule? 

Of course, one could ask a similar question about the probability of occurence of 
longer blocks of symbols after t iterations of the PCA rule. Due to the complexity of 
PCA dynamics, it is clear that questions of this type are rather hopeless if one wants 
to know the answer for an arbitrary rule. In spite of this, it may still be possible to 
provide the answer if the rule is sufficiently simple. 
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One of the methods which can be used to do this is studying the structure of 
preimages of short blocks and detecting patterns present in them. This approach has 
been successfully used by the author for a number of deterministic CA rules, such as 
elementary rules 172,142,130 (references 111 , 13 , and 13 respectively), and several 
others. It has also been used for a special class of PCA known as single-transition 
a-asynchronous rules Q. 

In this chapter, however, we would like to describe yet another method of com¬ 
puting probabilities of blocks of symbols, by setting up a system of recursive equa¬ 
tions which can then be explicitly solved. Such a recursive system can be easily 
constructed for any rule for probabilities of all blocks, but it is normally too big and 
to complex to be solved. In certain cases, however, one can find a smaller set of 
blocks for which the recursion is solvable. We will present one such example, using 
a PCA which can be viewed as a simple model of diffusion of innovations or spread 
of rumors. 

To give the reader a flavour of what to expect, let us informally define the afore¬ 
mentioned PCA rule. Suppose we have an infinite one-dimensional lattice where 
each site is occupied by an individual who has already adopted the innovation (1) or 
who has not adopted it yet (0). Initially the proportion of adopters is po. Once the 
individual adopts the innovation, he remains in state 1 forever. Individuals in state 
0 change can their states to 1 (adopt the innovation) with probabilities depending 
on the state of nearest neighbours: if only the right (resp., left) neighbour already 
adopted, the probability is p (resp., q), and if both of them already adopted, the 
probability is r. What is the proportion of adopters pt of after t iterations of the rule, 
assuming that the initial configuration is drawn from a Bernoulli distribution? We 
will show that the explicit formula for p, can derived. 


Pt = 


1 -E((poq-l){pop-l)y-F{l-ry if pqp^ - {p + q)po +r 
1 - {G + Ht){l-ry-^ if pqp^-{p + q)po + r = Q, 


where E,F,G,H arc constants depending on parameters p,q,r and po. 

In order to accomplish this, we will start from some general theoretical remarks, 
considering PCA as maps in the space of shift-invariant probability measures, sim¬ 
ilarly as done in Eoi, im, 0, and other works. More precisely, we will look at 
orbits of uniform Bernoulli measures under the action of PCA. 


2 Probabilistic cellular automata 

Probabilistic CA are often defined as stochastic dynamical systems. In this article, 
we will concentrate on Boolean CA in one dimension. Let s,(f) denote the state of 
the lattice site i at time f, where i G Z, f S N. We will further assume that i, (f) G 
{0,1} and we will say that the site i is occupied (empty) at time t if i,(f) = 1 (resp., 

Si{t)=0). 
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In a probabilistic cellular automaton, lattice sites simultaneously change states 
from 0 to 1 or from 1 to 0 with probabilities depending on states of local neigh¬ 
bours. A common method for defining PCA is to specify a set of local transition 
probabilities. For example, in order to define a nearest-neighbour PCA one has to 
specify the probability w(s,(f -b l))|s,-i (f) aKO Ai+i (0) that the site s,(f) with near¬ 
est neighbors (f) changes its state to s,(f -b 1) in a single time step. 

A more formal definition of nearest-neighbour PCA can be constructed as fol¬ 
lows. Let r be a positive integer, called radius of PCA, and let us consider a set of 
independent Boolean random variables A, b, where i G Z and b G {0,1}^''+^. Proba¬ 
bility that the random variable Xu, takes the value a G {0, 1} will be assumed to be 
independent of i and denoted by w(a|b). 


=a) = w(a|b). 


( 1 ) 


Obviously, w(l |b) -b w(0|b) = 1 for all b G {0,1}^''+'. The update rule for the PCA 
is then defined by 


Si{t + 1) 


( 2 ) 


Note that new random variables X are used at each time step f, that is, random 
variables X used at the time step t are independent of those used at previous time 
steps. 

Having the above definition in mind, we note that in order to fully define a 
nearest-neighbour PCA rule (i.e., rule with r = 1), it is enough to specify eight 
transition probabilities w{1\xiX2Xt,) for all xi ,X 2 ,X 3 G {0,1}. Remaining eight prob¬ 
abilities, w( 0 |xiX 2 X 3 ), can be obtained by w( 0 |xiX 2 X 3 ) = 1 — w(l |xiX 2 X 3 ). 

In any dynamical system, the main object of interest is the orbit of the system 
starting from a given initial point, and properties of this orbit. In the case of PCA, 
we often assume that the initial condition is “random” or “disordered”, typically 
meaning that each s, (0) is set to 1 with a given probability po and to 0 with prob¬ 
ability 1 — Po, independently of each other. We then want to answer question of 
the type “After t iterations, what is the proportion of sites in state 1 ?” or “After t 
iterations, what is the probability of finding a pair of adjacent zeros”? In order to 
pose and answer questions of this kind rigorously, we will present an alternative 
definition of PCA, as maps in the space of probability measures. 


2.1 Orbits of probability measures 

Let = {0,1} and X = . A finite sequence of elements of .e/, b = h\b2 

will be called a block (or word) of length ii. Set of all blocks of elements of jz/ of 
all possible lengths will be denoted by jz/*. 

A cylinder set generated by the block \i ~b\b 2 . ■ ■ ,b„ and anchored at i is defined 



as 


( 3 ) 
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The set of probability measures on the cr-algebra generated by cylinder sets of 
X will be denoted by S[ll(X). Details of construction of such measures, using Hahn- 
Kolmogorov theorem, can be found in 0. These details, however, are not essential 
for our subsequent considerations. Given a probability measure jJ. € ®T(X), measure 
of a cylinder set [b],, denoted by /r([b],), is often informally called a “probability of 
occurence of block b at site i”. 

Let the function w : x —[0,1], whose values are denoted by w(fl|b) for 

a G b G . 2 /^'’+*, satisfying ^e called local transition function 

of radius r, and let its values be called local transition probabilities. A probabilistic 
cellular automaton with local transition function w is a map F : 9ft(X) —S[ll(X) 
defined as 

(T^)([a],)= Y, w{a\h)p{[h]i-r) for alii G Z,a G.s/*, (4) 

where we define 

|a| 

w(a|b) = Y[w{aj\bjbj+i...bj+ 2 r). (5) 

j=i 

When the function w takes values in the set {0,1}, the corresponding cellular au¬ 
tomaton is called a deterministic CA. 

In this paper, we will exclusively deal with shift-invariant probability measures 
for which fi ( [b],) is independent of i. We will, therefore, drop the index i and simply 
write /4([b]). Moreover, we will be interested in orbits of Bernoulli measures 
defined for A € [0,1] by 

V;L([b]) (1forany be,e/*, (6) 

where #0(b) and #1 (b) denote the number of zeros and ones in b. In order to sim¬ 
plify the notation, we define 


P,(b) = (T'v;t)([b]), (7) 

which will be informally referred to as “probability of occurence of block b after t 
iterations of PCA rule F”. With this notation, eq. (0 can be written as 

Pt+iia)= Y vv(a|b)Pr(b), (8) 

for any a G si* and f S N. We will furthermore define 

Po(a) = V;i([a])=A''i(“)(l-Ar« (9) 


for any a G s/*. 

Elements of si* can be enumerated in lexicographical order, and corresponding 
probabilities arranged in an infinite column vector 
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P, = (P,(0),P,(1),P,(00),P,(01),P,(10),P,(11),P,(000)...)^. (10) 

Before we continue, note that not all these probabilities are independent. Due to 

additivity of measure, the following relationships, know as consistency conditions, 
are valid for any a € s/*, 

P,(a)=P,(a0)+B,(al)=P,(0a)+B,(la). (11) 

These conditions will be frequently used in our subsequent considerations. 

Since each P(+i(a), by the virtue of eq. (|8]l, is a linear combination of a finite 
number of Pr(b) values, we can write 

Pr+1=MP„ (12) 

where the infinite matrix M is defined by eq. (|8ll. This yields the following expres¬ 
sion for probabilities of all finite words. 


P, =M'Po, (13) 

where components of Pq are defined in eq. (|9]l. In theory, the above equation gives 
us a complete solution of the problem of determining the orbit of Bernoulli measure 
under iterations of a PCA rule. In practice, however, computing powers of an infinite 
matrix is a daunting, if not impossible, task. 

In practical applications, however, we rarely need all probabilities P;(a), that is, 
all components of the vector P,. Sometimes we are interested only in one specific 
probability, for example, T)(l). For a binary PCA, the expected value of a given 
lattice site after t iterations of the rule is equal to 1 ■ Pt{l) + 0 ■ Pt{Q) = Pt{l), and 
for that reason, P/(l) is sometimes referred to as an expected density of ones, to be 
denoted by pt, 

pt=m)- ( 14 ) 

Note that for Bernoulli measure Vx we have po = A. Given po, could one find an 
explicit expression for p, for a given PCA using eq. (fT2li ? This problem will be 
called a density response problem. Although it cannot be solved in a general case, 
we will demonstrate that for a sufficiently simple PCA it is a doable task. 

The idea is to setup a recursion similar to (fT^ . but using a “smaller” set of 
block probabilities, for which the matrix M has somewhat simpler structure, lending 
itself to direct computation of M'. If we could then express p, in terms of block 
probabilities from this “smaller” set, we would solve the density response problem. 

Let us define the concept of the “smaller” set first. A set of words s^* D 
C = {a 1,32,33,...} will be called complete with respect to a PCA rule F if for 
every a € C and f G N, Pr+i(a) can be expressed as a linear combination of 
P((ai),Pr(a 2 ),Pi(a 3 ),.... We will show a concrete example of a complete set in the 
next section. 


6 


Henryk Fuks 


3 Example PCA rule 

As an example, we will consider a PCA rule which generalizes some of the CA 
rules investigated in IT]. This PCA can be viewed as a simple model for diffusion 
of innovations, spread of rumors, or a similar process involving transport of infor¬ 
mation between neighbours. We consider an infinite one-dimensional lattice where 
each site is occupied by an individual who has already adopted the innovation (1) 
or who has not adopted it yet (0). Once the individual adopts the innovation, he 
remains in state 1 forever. Individuals in state 0 change can their states to 1 (adopt 
the innovation) with probabilities depending on the state of nearest neighbours. All 
changes of states take place simultaneously. This process can be formally described 
as a radius 1 binary PCA with the following transition probabilities, 

w(l|000) =0,w(l|001) = p, w(l|010) = l,w(l|011) = 1, (15) 

w(l|100) =^, w(l|101) = r, w(l|110) = l,w(l|lll) = 1, 

where p,q,r are fixed parameters of the model, p,q,r G [0,1]. In order to illustrate 
the difficulty of computing block probabilities for this rule, let us write eq. dH) for 
blocks a of length 1 and 2, 

p,+i(0) = P,(000) + (1 -p)T)(001) + (1 -^)Pr(lOO) + (1 - r)P,(101), 
P,+i(l)=pT-(001)+P,(010)+P,(011)+^P,(100)+rP,(101)+P,(110)-FP,(lll), 


P,+i(00) =P,(0000)-F (1(0001)-F (1(1000)-F (1-p)(l(1001), 
P,+1 (01) =pP,(0001)+ (l-p)P,(0010)+ (l-p)P,(0011)+p(l-^)P,(1001) 

-F (1 - r)P,(1010) -F (1 - r)P,(1011), 

P,+i (10) = (1 - ?)P,(0100) + (1 - r)Pj(OlOl) +^P,(1000) -F (1 - p)^P,(1001) 

+ (l-^)P,(1100) + (l-r)P,(1101), 

p,_^j(ll) = pP,(0010)-FpPr(0011)-F^Pr(0100) + rP,(0101)-FPr(0110)-FPr(0111) 
-Fp^Pr(1001)-FrP,(1010)-FrP,(1011)-F^Pr(1100) 

-F rPj(l 101)-FPKI 110)-FP,(1111). 

As we can see, in order to know Pt+i (1), we need to know probabilities of blocks of 
length 3 at time step f, and in order to compute these, we would need probabilities 
of blocks of length 5 at time step f — 2, etc. 

We will now show, however, that for the PCA rule defined in eq. (flSl l a complete 
subset of £/* can be constructed. This subset consists of clusters of zeros bounded 
by 1 on each side, that is, of blocks of the type 10"1, where n S N and 0" denotes n 
consecutive zeros. 

Proposition 1. The set C = { 101,1001,100001,...} is complete with respect to the 
PCA rule defined in in eq. ^751 ). 
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In order to prove this, we need to show that every ^r+l (10"1) can be expressed as 
a linear combination of probabilities of the type P,(10^1). Let us write eq. (l8]l for 
a = 10"1. Two cases must be distinguished, n= 1 and « > 1. For « = 1, we have 

P ;+1 (101) = p(l (01001)+ (1 -r)F;(01010) + (l -r)P;(01011) 

+ p^Fr(lOOOl) + (1 -p)qP,{l00l0) + (1 - p)^F;( 10011) + p(l -^)P;(11001) 

+ (l-r)Fr(l 1010)+ (l-r)Pr(l 1011). (16) 

By consistency conditions, P;(10010) + Pf(10011) = P;(1001) and f5(11010) + 
Pr(l 1011) =P,(1101), as well as F,(01001)+F,(11001) = P,(1001). This yields 

Ff+i(101) = (1 - r)F,(01010) + (1 - r)F,(01011) + pqPt{l000l) 

+ (l-p)^F,(1001)+p(l-^)F,(1001) + (l-r)F,(1101), (17) 

and further reduction is possible using F,(01010) + F((01011) = F((0101) and 
F,(0101)+F((1101)=F,(101).The final result is 

F;+i(101) = (l-r)F,(101) + (/7-2/7^ + ^)F,(1001)+MF,(10001). (18) 

For « > 1, using a similar procedure (omitted here), we obtain 

F,+i(10"l) = (l-p)(l-^)F,(10"l) + (p-2/7^ + ^)F,(10"+*l)+p^F,(10"+2l). 

(19) 

Equations (fTsT l and (fTOl l clearly show that the set C is complete. □ 


4 Calculations of (10" 1) 

Having a complete set of block probabilities, we can now write eq. dUll and dull in 
matrix form. 


1 

li 

1_ 


aficOOO--- 
0 a b c 0 0 


'Pt(lOl) ' 
Pi(lOOl) 


_ 

0 0 a b c 0 



7’*+i(10"l) 


0 0 0 a b c 


Fj(10"l) 







where a = {1 — p){l — q),d = 1 — r, b = p — 2pq + q, and c = pq. 
Let us define 
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d b c 0 0 0 ■■■ 


■Pr(lOl) ■ 

0 a b c 0 0 


Pr(lOOl) 

0 0 a b c 0 

, Pr = 


0 00 a b c 


P,(10"l) 

- ■ 




( 21 ) 


We will use diag(xi,X 2 ,X 3 ,...) to denote an infinite matrix with xi,X 2 ,X 2 ,,. ■ ■ on 
the diagonal and zeros elsewhere. Similarly, sdiag( 2 Ci,X 2 ,X 3 ,...) will denote shifted 
diagonal matrix having xi,X 2 ,X 3 ,... on the line above the diagonal and zeros else¬ 
where, and ^sdiag(xi ,X 2 ,X 3 ,...) will denote doubly-shifted diagonal matrix, with 
xi ,X 2 ,X 3 ,... on the second line above the diagonal and zeros elsewhere. With this 
notation, we have 


where 


Now, 


M = A + B + C, 

(22) 

A = diag(5,a,a,...), 

(23) 

B = sdiag{b,b,b,...), 

(24) 

C = ^sdiag(c,c,c,...). 

(25) 

P,=M'Po, 

(26) 


and we need to compute M'. We will do it by considering a special case first. 


4.1 Special case: a = a 


When a = a, matrices A, B, and C pairwise commute, thus we can use the trinomial 
expansion formula, 

M' = (A + B-FC)'= Y. ( ■ ^ (27) 


where 

Generalizing the previously introduced notation, let 


(28) 


"sdiag(xi,X2,X3,...) 


(29) 


denote n-times shifted diagonal matrix, which has xi,X 2 ,X 3 ,... entries on the n-th 
line above the diagonal and zeros elsewhere. It is straightforward to prove that 
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A'= diag(a',a',a',. ■ ■)) (30) 

B' = hA\d,g{b\b^, (31) 


and, consequently, 

A'B^C'^ =2+2^ sdiag(flW,a'Vc^flW,.. .)■ (33) 

In the first row of the above matrix, the only non-zero element (a'b^c'^) is in the 
column 1 -I- y -I- 2k. In the second row, the only non-zero element {a'b^(^) is in the 
column 2 -I- y -I- 2k, and so on. This means that 


A'B^'c'^Po = 


a'l72c^fb(10^+2+2^1) 
a'b^ c'^Pq{\(P-+}+^H) 
a'y(^Po\\(f’+^+^h) 


(34) 


Using the above and the fact that Po(10"l) = Po (1 — Po)”, we can now write 

■d'VcVo(l-po) 1+2+2^ 
a'VcVo(l -Po) 2 + 2 + 2 ^ 

a'VcVo(l 


P, = M'P„= £ 

i+i+k=t 


We finally obtain 

Pt(io'i)= Y. 


i+i+k=l 

= Po(l-Po)' L 


«VcVo'(l-po)'+^+''^ 

t 


a‘[b{l-po)Y[c{l-poy] 


2-ik 


(35) 


i+j+k=t 

= p^(l -po)' (fl-bfe(l-po)+c(l -Po)^)'. (36) 


4.2 General case 


We are now ready to handle the general case, without the a = a assumption. Let us 
first note that f-th powers of matrices 
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d b c 0 0 0 ■■■ 

t 

a b c 0 0 0 ■ ■ ■ 

0 a b c 0 0 


0 a b c 00 

0 0 a b c 0 


0 0 a b c 0 

0 0 0 a b c 

5 

0 0 0 a b c 





differ only in their first row. This implies that the expression for P, (10^ 1) given in eq. 
(|3^ remains valid for I > 1 even if a ^ a. We only need to consider / = 1 case, that 
is, to compute Pr(lOl). This can be done by writing eq. (fTSl) and replacing P,(1001) 
and P,( 10001) by appropriate expressions obtained from eq. (l36]) . 

P ,+1 (101) = aP,{m\)+bpl{\- Po)" {a + b{\ - Po) +c(l - po)")' 

+ cp^(l -po)^ (a + fe(l-po)+c(l -po)^)'. (38) 


This can be written as 


P,+ i(lOl)=flP,(lOl)+/:0', (39) 

where 

^ = ^Po(l-Po)^ + cPo(l-Po)^ (40) 

0=a + /7(l-po)+c(l-po)^ (41) 

Eq. (l39l l is a first-order non-homogeneous difference equation for P((101), and, as 
such, it can be easily solved by standard methods El. The solution is 

t 

P,(101) =Po(lOl)5'+/:^d'^'0'^‘. (42) 

!=1 

The sum on the right hand side is a partial sum of geometric series if 5 7 ^ 0, or of 
an arithmetic series when 5 = 0. Using appropriate formulae for partial sums of 
geometric and arithmetic series one obtains 


Pr(lOl) 


Po(ioi)5'+/:(fl'-0')/(5-0) ifd^e, 
Po{lQl)a^ + Kd‘^^t ifd = 6. 


(43) 


Taking Po(lOl) = p^(l —p) and replacing K and 0 by their definitions we obtain, 
for a — a + b{l — po) +c(l — po)^ ^ 0 (which is equivalent to 5 7 ^ 0 ), 


fl-a + b(l — Po)+c(l — po)"^ 

^ Po^(l-po)(a-a) ^ 

a-d + b{l -po) +c(l -po)^ 


( 44 ) 


Pr(lOl) 
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For a — a + b{l — po) + c{l — po)^ = 0 (equivalent to a = 0), the solution is slightly 
simpler, 

P;(101) = Pq (1 — Po) ((cpo^ — (fe + 2c)Po + + c)f + fl)5^ *■ (45) 

We now have expressions for P,(10^1) for / = 1 (ea.l44landl45]) and for I > 1 (eq. 

[36ll. 


5 Cluster expansion 


We are finally ready to compute Pi. To do this, we will use the formula 

P,(0) = f;^F,(10'^l), (46) 

1=1 

which we will refer to as “cluster expansion”. Various proofs of this formula can be 
given (see, for example, ifT^ l. but we will show here that it is a direct consequence 
of additivity of measure. 

Consider a cylinder set of a single zero anchored at i, [0],. A single zero must 
belong to a cluster of zeros of size k with possible values of k varying from 1 to 
infinity. If it belongs to a cluster of k zeros, than it must be the y-th zero of the 
cluster, with possible values of j varying from 1 to k. Therefore, 

[0L=UU[10'W (47) 

k=\i=\ 

Since all the cylinder sets on the right hand side are mutually disjoint, their measures 
add up, thus 

oo k 

(7^'v;t)([O]0 = E E(^’'v;i)([10'l],•-;•)• (48) 

k=\i=\ 

The measure is shift-invariant, thus [10^1 ],_;■) = F?(10^1), and we obtain 

7^r(0) = EE7^/(10'l)> (49) 

k=\i=\ 


which yields eq. (|46] |, as desired. 

We can now compute PtiQ) using the cluster expansion formula and eq. (l44l i. (l45l l 
and (l36l l. We will first consider the case of a — a + b{l — po) -|-c(l — po)^ ^ 0, that 
is, using eq. (l44l i forP,(101). 
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p,i0) = tmo'l) = {a + bil-po)+cil-pofy 

a-a + b{l-pQ)+c[l-poy ^ 

pO^{l-po){a-a) ^ ^ ^ ^ 2 ^^ 




Since 


we obtain 


;ri ^ \i (l+Po)(l Po)^ 

—’ 


(50) 

(51) 


P,(0) = 


Po (1-po) (P + c-cpo) 

fl-a + P(l -po) + c(l -po)^ 


(a + b{l — Po) +c(l — Po)^) 


+- .. «' + (1 +Po)(l -Po)' (« + ^7(1 - Po) + c(l - Pof )’. 


!-fl + /7(l -Po)+c(l -po)2 


(52) 


After substitution of d,a,b,c and simplification, as well as taking Pr = 1 —Pt{Q), the 
following expression for p, is obtained. 


pqpo -{p + q)Po + r 

pqpo^ - (p + q)po + r 

When a — d + b{l — po) + c( 1 — po)^ = 0, similar calculations can be performed, 
but this time using using eq. ( l45l l for P((101)- After simplification, this yields 


p, = 1 - (1 - Po) (po (po - 1) (pqPo -p + pq-q)t + l-r){l-r) 


r-l 


(54) 


Let us summarize this in a more readable form, noticing that after substitution of 
d,a,c,bhy their definitions the condition a — fl+P(l — po) + c(l — po)^=0 becomes 
PqPo — {p + q)Po + f — 0. Our final expression for the density of ones can be written 
as 


Pt = 


l-E{{pQq-l){pop-l)y-F{l-ry if pqp^ - {p + q)po +r ^ 0, 

1 - (G + //f)(l -r)'^* if p^p^- (p + ^)po + r = 0, 

(55) 


where definitions of E,F,G,H can be figured out by comparing the above to eq. (l52] | 
and (l53]) . 

We can see that in the non-degenerate case (when pqp^ — {p + q)po + r 7^0), 
the limit poo = lim,^ooPr always exists, and that p, approaches p„c exponentially 
fast, excluding special cases when p, = const (such as po = 0 or p = o' = 0, r = 1). 
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In the degenerate case, always exists as well, but the approach of pt to is 
linearly-exponential. 

It is worth noting that the existence of the degenerate case is a fairly subtle phe¬ 
nomenon, and that it would be very difficult to discover the linearly-exponential 
convergence by computer simulations alone. This illustrates the point that having 
a formula for p, brings some advantages, and that the search for such formulae is 
worthwhile. 

As a separate remark, let us note that deterministic CA are nothing else but spe¬ 
cial cases of PCA, thus we can choose integer values of p,q,r and obtain relevant 
expression for pt for a number of elementary CA rules (ECA), as follows. 

ECA rule 206 (p = l,q = 0,r = 0) or rule 220 (p = 0,q= l,r = 0) 

p, = l-po(l-po)-(l-po)'+^ (56) 


ECA rule 222 (p = q = l,r = 0) 

„ , , 2(l-po)2'+'-fpo(l-po) 

Po —2 

ECA rule 236 (p = 0,^ = 0,r = 1) 

p, = l-(po^-l)(l-po)^ 

ECA rule 238 {p = l,q = Q,r = 1) or rule 252 (p = Q,q = l,r = 1), 

p, = l-(l-po)*+i, 


(57) 


(58) 


(59) 


ECA rule 254(p = ^ = r= 1) 


p, = l-(l-po)2'+^ 


The above formulae agree with those derived informally in JT] . 


(60) 


6 Conclusions 

We presented a method for computing the density of ones in the orbit of the 
Bernoulli measure under the action of a probabilistic cellular automaton, using a 
simple PCA rule as an example. Eor this rule, we were able to construct a complete 
set of block probabilities, and then solve the resulting recurrence relations. By using 
the cluster expansion, we then obtained the required density of ones. Although this 
method is obviously applicable only to PCA rules with rather simple dynamics, it 
may be possible to find other PCA rules with complete sets, thus making the method 
useful for them. Generalization of the rule used in this paper to larger neighbour¬ 
hood sizes comes to mind as a first possibility, and sufficiently simple deterministic 



14 


Henryk Fuks 


rules, such as asymptotic emulators of identity investigated in ID, are also possible 
candidates. 
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